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M.  ASCE 


INTRODUCTION 


In  this  paper,  the  minimax  method  is  applied  to  boundary  value  problems 
which  arise  in  structural  mechanics.  This  is  a  weighted  residual  method  in 
which  the  maximum  absolute  value  of  a  residual  is  minimized.  Like  other 
weighted  residual  methods  (7,  13),  a  trial  function  is  employed  which  con¬ 
sists  of  undetermined  parameters  and  basis  functions.  This  trial  function 
is  introduced  into  governing  differential  equations, and  the  maximum  absolute 
residual  among  several  residuals  at  discrete  points  in  the  domain  is  minimized. 
This  residual  minimization  criterion  is  applied  to  a  finite  element  formula¬ 
tion  by  using  piecewise  trial  functions  defined  on  each  element. 

Computational  implementation  is  achieved  using  linear  programming.  This 
approach  has  been  used  in  the  global  sense  to  obtain  solutions  for  differen-  1  ' 
tial  equations  (3,  13,  21,  28). 

Since  the  linear  programming  technique  can  be  applied  to  solve  over- 
determined  systems  of  equations  (24),  more  mesh  points  than  the  number  of 
unknown  parameters  can  be  used  to  improve  the  solution.  This  is  in  contrast 
to  the  collocation  method.  Also,  equality  constraints  representing  the 
boundary  conditions  and  inter-element  continuity  conditions  can  be  included 
in  the  formulation.  This  feature  provides  more  freedom  in  choosing  a  trial 
function,  which  is  often  the  most  important  step  in  the  weighted  residual  method. 
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In  the  collocation  method,  the  roots  of  the  orthogonal  polynomial 
are  frequently  used  as  collocation  points  with  considerable  success  (10, 

23,26,27).  This  is  referred  to  as  the  orthogonal  collocation  method.  Convergence 
studies  (3,10,12,20)  of  the  orthogonal  collocation  method  show  that  in  most 
cases  a  convergence  rate  similar  to  that  for  the  Galerkin  and  Ritz  methods 
can  be  obtained.  This  procedure  will  be  utilized  here  in  choosing  optimal 
locations  for  mesh  points. 

Since  inequality  constraints  are  acceptable  in  a  linear  programming 
problem,  the  minimax  method  can  be  used  for  limit  analyses  and  can  analyse 
structures  with  off-set  supports,  which  pose  a  contact  problem.  When  multiple 
optimal  feasible  solutions  arise  for  these  types  of  problems,  a  parametric 
programming  method  described  in  Appendix  I  can  be  used  to  choose  a  specific 
optimal  solution. 

Soon  after  the  upper  bound  theorem  and  the  lower  bound  theorem  in 
plasticity  became  available  (19)  and  the  linear  programming  technique  was 
developed  by  Dantzig,  limit  analysis  was  identified  as  a  linear  programming 
problem-  and  linear  programming  algorithms  have  been  used  for  limit  analysis 
solutions  since  then  (1,4  ,5,10,15,17,29,30). 

This  method  appears  to  be  easy  to  formulate  and  simple  to  use.  Like 
the  collocation  method,  the  proposed  method  does  not  involve  the  integra¬ 
tion  that  is  necessary  with  the  other  weighted  residual  methods. 


FORMULATION 

Elasto-Static  Analyses.-  For  a  given  governing  differential  operator 
equation  of  motion 

Au  =  f  (1) 

we  seek  an  approximate  solution  u  whose  image  under  differential  operator 
A  approximates  the  given  function  f.  Let 


I 

u(x)  =  X  a  (x) 
-  i-1  1  1  " 


(2) 


i=l 
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In  the  proposed  method,  the  criterion  to  determine  the  unknown  para¬ 
meters  a.  is  to  minimize  the  maximum  absolute  residual  among  several 
residuals  at  discrete  points.  Let 

r  =  max |  R(Xj ) !  j  =  1,2,.  .  .  .  J  (4) 

in  which  x  are  the  locations  of  mesh  points.  Then,  a  linear  programming 
problem  can  be  established  (24).  Find  the  parameters  a.  such  that  Z  =  r 
is  minimized  and  the  constraints 


R(x.)  -  r  <  0 

3  j  =  1,2,.  .  .  J 

-R(Xj)  -  r  <  0 

are  satisfied. 


(5) 


If  the  trial  function  does  not  satisfy  some  of  the  boundary  conditions, 
the  following  equality  constraints  can  be  included  in  the  above  formulation 
(Eg.  5). 


Bu(x  )  =  b(x  )  p  =  1,2,.  .  .  .  P  (6) 

-P  ~P 

in  which  B  is  a  boundary  differential  expression,  b(x)  is  a  given  function 
defined  on  the  boundary,  and  x^  are  the  locations  of  boundary  mesh  points. 

When  the  finite  element  method  is  used,  interelement  continuity 
conditions  should  be  satisfied.  That  is,  if  p  is  the  number  of  the  highest 
order  appearing  in  the  governing  differential  equation  ,  the  approximate 
function  and  up,  to  (p-1)  derivatives  must  be  continuous  across  interelement 
boundaries  (31).  In  the  minimax  method,  these  conditions  can  be  placed 
as  constraints  in  the  linear  programming  formulation.  However,  the  high 
order  inter-element  continuity  requirement  tends  to  increase  the  problem 
size  by  increasing  the  number  of  nodal  variables  or  by  usinn  additional 
constraint  equations.  This  constitutes  a  drawback  of  the  proposed  mini- 
max  method  compared  to  conventional  finite  element  methods. 

If  the  number  of  constraints  is  greater  than  or  equal  to  the  number 
of  unknown  parameters,  the  existence  and  uniqueness  of  a  polynomial  of 
best  approximation  is  guaranteed  (9,  22). 

Limit  Analysis.-  According  to  the  lower  bound  theorem  in  plasticity, 
if  an  equilibrium  distribution  of  stress  can  be  found  which  balances  the 
applied  load  and  is  everywhere  below  yield  or  at  yield,  the  structure 
will  not  collapse  or  will  be  just  at  the  point  of  collapse  (19). 

To  apply  the  previously  formulated  minimax  method  to  limit  analysis 
by  using  the  lower  bound  theorem,  the  trial  function  should  satisfy  the 
equilibrium  equations  in  the  domain  and  boundary  conditions  on  the  boundary, 
pattern  of  load  distribution  is  known.  However,  the  magnitude  of  the  load 
factor  A  is  unknown  and  is  treated  as  an  unknown  variable  together  with 
unknown  coefficients  of  a  trial  function.  The  yield  condition  is  introduced 


The 


at  some  check  points  as  inequality  constraints.  The  minimization  of  the 
maximum  residual  produces  only  an  equilibrium  state  for  any  value  of  A. 

In  other  words,  for  any  A  which  does  not  violate  the  yield  condition,  the 
values  of  unknown  coefficients  can  be  determined  for  the  resulting  stress 
distribution  to  balance  the  applied  load,  A,  accordingly.  This  means  that 
the  linear  programming  problem  has  multiple  solutions  unless 
there  is  another  constraint  that  can  specify  A.  Our  purpose  is  to  obtain 
a  solution  that  has  the  maximum  value  for  A  among  these  multiple  solutions . 
Therefore,  we  have  to  minimize  the  maximum  residual  and,  at  the  same  time, 
maximize  A.  Then,  the  linear  programming  becomes:  find  the  a.  such  that 
2  =  r  is  minimized,  \  is  maximized,  and  the  constraints 


R(Xj)  -  r  <_  0 
-R(x.)  -  r  <  0 
Bu(x  )  =  b(x  ) 

-p  -p 

Y  (x  )  <  k 
—n  — 


j  1,2,.  .  .  J 


p  =  1,2, .  .  .  P 
n  =  1,2,  .  .  .  N 


are  satisfied.  Here  Y(xn)  is  the  value  of  the  yield  function  evaluated 


at  x  =  x  . 
—  — n 


This  problem  of  multiple  objective  functions  can  be  handled  as  a  linear 
programming  problem  using  the  parametric  programming  method  described  in  Appendix 
I.  The  technique  in  Appendix  I  of  obtaining  a  solution  with  a  maximum  A  is  to 
use  a  new  objective  function  Z  =  r  -  eA,  where  £  is  a  small  positive  number. 

The  solution  gives  a  lower  bound  to  the  limit  load. 


NUMERICAL  EXAMPLES  AND  RESULTS 

Beam  on  Elastic  Foundation.-  Consider  a  beam  on  an  elastic  foundation 
subject  to  a  uniform  load  as  shown  in  Fig.  1. 

The  governing  differential  equation  for  engineering  beam  theory  is 

EIy1V  +  kgy  =  q  (3) 

in  which  E  is  Young's  modulus  of  elasticity,  I  is  the  moment  of  inertia 
of  the  beam  cross  section,  ks  is  the  stiffness  of  the  elastic  foundation, 
and  q  is  the  uniform  load  intensity. 

First,  an  ordinary  polynomial  is  used  as  a  trial  function. 

I  .  , 


y(x)  =  Z  a.  x 
,  ,  1 
i=l 


Then, 


R(x)  =  Ely  (x)  +  kgy(x)  -  q  (10) 

Let  one  element  apply  for  the  entire  beam.  Use  an  equidistant  spacing 
mesh  for  the  numerical  calculations.  The  results  shown  in  Table  1  are 
quite  accurate. 
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An  exception  occurs  for  the  beam  of  length  L  =  2  with  a  10th  order 
polynomial.  In  the  use  of  a  very  high  order  polynomial,  the  differences 
between  magnitudes  of  the  coefficients  in  the  constraint  equations  become 
very  big  and  significant  numerical  error  results.  This  is  a  drawback  in 
using  an  ordinary  polynomial.  However,  this  difficulty  can  be  overcome  by 
refining  elements.  An  advantage  of  using  an  ordinary  polynomial  is  that 
there  is  no  need  for  matrix  inversion  to  transform  the  coefficients  of  the 
polynomial  to  nodal  variables  as  is  encountered  in  the  usual  finite  element 
formulations  (31). 

n  Next  use  Hermitian  shape  functions  as  bases  of  the  trial  function.  Let 
H. .  (£)  denote  the  nth  order  Hermitian  shape  function  which  is  the  (2n+l)th 
or3er  polynomial  (31).  In  this  polynomial,  n  is  the  number  of  derivatives 
that  the  set  can  interpolate,  i  is  the  order  of  derivatives  of  II.  .  (Q  with 
respect  to.,  £,  and  j  =^1  or  2  are  the  element  node  numbers.  In  numerical 
example,  Ht .  (Q  and  H..(^)  are  used.  Then,  the  trial  function  for  each 

r  13  13 
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element  becomes 
2 

9„ 


or 


ih  «*<«>»,  ♦  <j«  r,B*4i'3  *.5a> 


y^  ■  I  (fr.(;)  y.  +  Hl\(c)  y  +  H~.(H  y  __  +  H^.(0  y  _  )  (1 

e  4.i  03  ID  *£)  23  *55j  33  'S^j 


(11) 


2) 


respectively. 

As  shown  in  Table  2,  the  numerical  results  are  quite  accurate  and, 
for  the  same  number  of  unknown  >,  increasing  the  number  of  mesh  points 
improves  the  solution.  It  is  also  noted  that  improvement  can  be 
achieved  by  refining  elements. 


When  there  is  a  concentrated  load,  the  structure  can  be  divided  into 
finite  elements  such  that  the  location  of  the  concentrated  load  becomes 
a  node  point.  The  discontinuity  of  a  variable  such  as  y"’  should  be 
taken  into  account  in  the  formulation.  This  technique  is  also  used  in 
the  analysis  of  beams  with  off-set  supports. 

Torsion  of  Prismatic  Bar.  -  As  an  application  of  the  minimax  method 
to  partial  differential  equations  consider  the  elastic  torsion  of  a 
prismatic  bar.  The  governing  differential  equation  (25)  is 


in  the  domain 


(13) 


0=0  on  the  boundary 


(14) 


in  which  <t>  is  Prandtl's  stress  function,  G  is  the  shear  modulus,  and 
6  is  the  angle  of  twist  per  unit  length  of  the  bar. 

Since  the  highest  order  derivative  appearing  in  the  governing 
differential  equation  is  two,  the  requirement  of  interelement  continuity 
is  the  continuity  of  0,  $,x,  «J>,  .and  <p,  .  Therefore,  the  following  trial 

function  can  be  used  for  tne  rectangular  element  shown  in  Fig.  2, 


V?'nl  •  ,fx  I»oi<5)  »oV"’  *ij  *  hVs’  “oV™  ».?il 

*hJ.U)  Hj.(n)  *,nij  ♦  »;.,?)  H^ln)  ♦,e|i  1 


The  residual  in  each  element  is 

1  ^e  1  d\ 

V-3'n)  =  I2  TP~  +  b2  +  2Ge 


(15) 
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As  a  numerical  example,  a  prismatic  bar  of  a  square  cross  section 
will  be  treated  by  the  minimax  method.  Let  H  be  the  height  and  W  the 
width  of  the  cross  section.  The  analytic  solution  for  the  torsion  of  a 
rectangular  bar  as  taken  from  Ref.  75  is 

T  =  KG(?W 
max 

M  «  K1Gew3H 

in  which  T  is  shear  stress,  is  the  torque,  and  K  and  are  numerical 
factors  depending  on  the  ratio  H/W. 

Due  to  the  symmetry,  only  one  quadrant  of  the  cross  section  needs 
to  be  considered,  results  are  shown  in  Table  3.  For  the  same  number  of 
mesh  points,  the  selection  of  mesh  points  at  the  Gaussian  points  gives 
better  results  than  the  use  of  mesh  points  at  a/4  and  2a/4.  It  is  also 
noted  that  the  solution  for  K  with  3x3  Gaussian  points  is  better  than 
with  2x2  Gaussian  points.  However,  2x2  Gaussian  points  plus  a  center 
mesh  point  gives  poorer  results  than  simple  2x2  Gaussian  mesh  points. 
This  indicates  that  convergence  with  respect  to  increasing  mesh  points 
is  not  monotonic.  The  results  also  indicate  that  increasing  the  number 
of  elements  improves  the  results  even  with  a  modicore  choice  of  mesh 
points . 

Plane  Stress  Analysis.-  For  a  plane  (x,y)  elasticity  problem  the 
equilibrium  equation  is 


3a 


XX 


3x 


?o 

XV 

3y 


*  0 


3a 

-JZ  + 

o  y 


3o 


3x 


^  0 


(16) 


if  body  forces  and  thermal  effects  are  neglected.  In  terms  of  displacement 
components  u  and  v  ,  this  becomes 


1  1+u  1 

1-v  U'xx  2(l-v)  v'xy  2  U'yy 


0 


(17) 


1  1+v  1 

1-v^'yy  2 ( 1- v)  u'xy  2  v'xx 


0 


The  solution  to  Eqs  (17)  subject  to  appropriate  boundary  conditions  consti¬ 
tutes  the  solution  of  the  problem  of  elasticity.  Since  the  highest  order 
of  derivatives  appearing  in  the  governing  differential  equations  is  two  in 

Eq .  (18),  u,  u,  ,  u,  ,  and  u,  should  be  continuous  between  elements.  The 
x  y  x> 


sane  applies  to  the  displacement  component  v.  Thus,  a  conforming  element 
using  the  Hermitian  shape  functions  as  basis  functions  is  used. 

'“I 

u  =  I  I  [H1  C ;,)  H*  (p)  u.  .  +hJ-.(.)  H*.u,  . 

i=l  j  =  l  01  °3  11  °3  cl3 

+  hJ  .(0  h  J  (n)  u, ,  .  +  h]"  .  (^)  hJ-.(h)  u,  . 

03  i]  "13  ii  s  13  '  013 


2  2 


~  i  [hJ;.(-:)  hJ-.(i)  v..  +h.1.(-)  hJ.h)  v _ 

i=1  j=1  01  °i  ij  ix  j  03  a: 


+  HJiU)  Hij(n)  v'  nij  +  Hii(C)  Hij (n)  v'-nii 


in  which  u  ,  v  are  trial  functions  for  u  and  v,  respectively,  in 


element . 


As  an  example  of  a  plane  stress  problem,  consider  a  plate  subject 
to  simple  in-plane  forces  such  as  pure  tension,  pure  bending,  and  pure 
shear.  For  the  pure  tension  problem  shown  in  Fig.  3,  the  approximate 
boundary  conditions  are 

* 

U1  =  U2  =  U3  =  11 '  yl  =  U'y2  =  U'y3  =  U'y4  =  u'y7 

=  V!  =  v4  "  V7  =  V'xl  =  V ' x2  =  V'x3  =  V ' x4  =  v'x7  =  0 


j  =c  =0  =  m 

xx7  xx8  xx9 


yy3  yy6  _  xy6  °xy8  xy9 

in  which  u^  is  the  value  of  u  at  node  1,  u,  is  the  value  of  — ^  at 
node  1,0  _  is  the  value  of  O  at  node  7, "etc.  °*/ 

XX  /  XX 

The  results  obtained  with  the  minimax  method  for  the  simple  in-plane  forces 
are  exact  and  this  formulation  passes  the  so-called  patch  test.  According 
to  this  test,  in  order  for  a  solution  to  converge  to  the  correct  one  by 
refining  elements,  a  patch  of  elements  subjected  to  a  specific  nodal  dis¬ 
placement  corresponding  to  a  state  of  constant  strain  should  produce  the 
constant  strain  state  throughout  the  elements  (6) . 

Limit  Analysis  of  a  Fixed-fixed  Beam.-  Suppose  the  limit  load  is  sought 
for  a  prismatic  fixed-fixed  beam  subject  to  a  uniform  load  as  shown  in 
Fig.  4.  Due  to  the  symmetry,  only  half  of  the  beam  needs  to  be  considered. 
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Let 


y(^)  =  I  [Hoi(^>  yi  +  H  (Q  y  r.  +  H  (0  y  +  H  At,)  y  rr 

A  AJ*  #  v»i  -«*■  f  t>S>i  ^  r>  -»•*■ 

For  a  beam  problem,  the  yield  condition  is  expressed  as 


rrrj  (19, 


|M  I  <  M  (20) 

1  n 1  —  p 

in  which  is  the  moment  at  x  =  xr  and  M  is  the  plastic  yield  moment 
of  the  beam  which  is  obtained  when  the  whole  section  of  the  beam  becomes 
plastic  (2),  as  shown  in  case  C  in  Fig.  5.  If  the  yield  condition  is 
checked  at  £  =  0,  5=1 


M 

n 


EX 

I7  y'55n 


n  =  1,  2 


(21) 


Choose  the  uniformly  distributed  load  q  =  1.  Then,  the  residual  is 


R(5) 


d5“  El 


(22) 


The  linear  programming  problem  becomes: 
that  2  =  r  is  minimized,  \  is  maxmimized 


Find  y.,  y,,.,  y,^., 
,  and  the  constraints 


y 


. such 


R(-t  )  -  r  <  0 

- 

R(^j)  -  r  <_  0 


El 


<  M 


c 


<  M 

-  p 


(23) 


are  satisfied.  Using  g.  =  0.2,  0.4,  0.6,  0.8,  and  1,  the  following  results 
were  obtained  J 


Ml  =  -12  lb-in  (1356.36  N-mm) ,  =  7.2  lb  (32.04  N) ,  y2  =  0.004  in  (0,1016  mm) 

Mn  =  6  lb-in  (678.18  N-mm),  A  =  1.44  lb  (6.408  N) 


Here , 


M.  and  V.  are  the  moments  and  shear  force  at  node  i. 

l  l 


1.44  is  the 


10 


lower  bound  on  the  limit  load  factor  and  causes  the  plastic  moment  at 
=  0. 

If  the  material  is  elastic  perfectly  plastic,  remains  the  same 
as  M  .  Upon  a  further  increase  in  load,  the  elastic  part  of  the  beam 
will^supcort  the  increase  in  load.  Therefore,  by  solving  the  following 
linear  programming  problem,  the  next  largest  load  which  causes  further 
yielding  can  be  determined:  Find  v  ,  y>r.,  y,  rQ»  and  y,  rrr.  such  that 
Z  *  r  is  minimized,  X  is  maximized,  and^lhe  constraints 


R(  ^ . )  -  r  <  0 
:  - 

-R(  L)  -  r  £  0 

(24) 

EI 

-  p-j  y,  =  Mp 

are  satisfied.  The  numerical  results  for  this  problem  were  found  to  be 

M1  =  -12  lb-in  (-1356.36  N-mrn)  =  9.6  lb  (42.72  N) ,  yn  =  0.11  in.  (2.794  mm) 
M„  =  12  lb-in  (1356.36  N-mm)  X  =  1.92  lb  (8.544  N) 

As  indicated  in  Fig.  4(d),  the  collapse  mechanism  has  been  formed  and 
\  -  1.92  is  the  exact  load  factor  for  the  collapse  load. 


CONCLUSION 

It  is  shown  that  the  minimax  weighted  residual  method  can  be  used  for 
obtaining  finite  element  solutions.  The  method  appears  to  be  relatively 
easy  to  set  up  and  gives  satisfactory  results  for  the  example  problems. 

This  method  is  very  attractive,  particularly  for  ordinary  differential 
equations  and  low  order  partial  differential  equations.  It  can  be  used  to 
solve  problems  with  inequality  constraints. 

Unlike  the  collocation  method,  the  solution  can  be  improved  by  using 
more  mesh  points  than  the  number  of  unknown  coefficients  in  the  given  trial 
function.  Also,  the  solution  can  be  improved  by  refining  elements.  The 
Gaussian  points  are  optimal  mesh  points  for  the  proposed  minimax  method. 
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PARTICULAR  SOLUTION 


APPENDIX  I  -  PARAMETRIC  P ROGRAM-MING  METHOD  TO  CHOOSE  A 
FROM  MULTIPLE  SOLUTIONS 


f  T 

Consider  the  problem,  find  'xl  such  that  Z  =  c!  {x/  are  minimized 
and  the  constraints 

[A ] {x}  =  {b} 

i'xf  >  0 

are  satisfied.  Let  some  perturbation  be  given  to  c 


rC+}T  =  {c}T-  £{f I' 


in  which  £  is  small  positive  number.  Mow,  assume  that  {x^}  is  one  of  the 
multiple  solutions  {x^}.  Then,  2+  becomes  (14) 

’-+  ■  -  «<«?’>  <v'  -  fe?T  <*»  -  *c^r 

If  there  are  l  multiple  solutions 


{cB}  (XB}  =  {CB}  {X3}  = 


T 

Qi  r  Ci 


=  fcj)  {*?> 


'C  b  ‘  lV 


If  we  can  have  one  of  x  remain  optimal,  minimizing  Z  means  also  maximizing 
,  T 

c{f0-'  {x  and  tnis  is  the  same  as 


*  r  T 

Maximize  Z  =  c{f}  {x> 

such  that  [ A ] { x }  =  {b } 

{x}  £  lx  } 

B 


Therefore,  by  minimizing  Z  =  {c  }T{x}  such  that  [A]{x}  =  {b},  {x}  _>  0 
and  by  properly  choosing  f.,  we  can  obtain  a  particular  solution  from  the 
multiple  solutions  which  satisfy  min  Z  =  {c},;,{x}  and  max  Z*  =  £{f  }T{x}  such 

that  [a]{x}  =  {b},  {x^6{Xg},  e  can  be  chosen  as  follows: 

We  wish  to  maintain  one  of  the  {x0}  as  the  optimal  basic  solution  for  the 
new  problem 


{c+ }T {x  > 


minimize  Z  = 
such  that  [A] (x)  =  (b } 

{ x }  >  0 


denote  by  zt  -  ct  the  value  at  z . 

j  3  J 


Cj  when  {c }  is  replaced  by  {c  }  , 


then  the  critical  value  of  e  is  such  that  any  increase  in  £  would  make  one 
+  + 

or  more  -  c^  positive 


2j  ~  Cj  =  ({CB}T  “  £{fB}){yj}  *  Cj  "  £fj  =  2j  “  Cj  -  ' ( { f B } 


“  f,) 


in  which  fg  is  the  row  vector  that  contains  the  components  of  f  corresponding 

T 

to  the  components  of  c  in  cg.  If  {fg}  {y^ }  -  f ^  are  nonnegative,  then  we 

can  make  £  arbitrarily  large  without  destroying  optimality.  However,  if 
T 

one  or  more  {.f  }  (y .  }  -  f.  are  negative  and  e  is  large  enough,  the  corresponding 
B  j  1 

+  +  J 

z.  -  c.  will  become  Dositive .  Thus,  the  critical  value  of  £  >  0  is  given  by 

:  3  c  - 


or 


min 


z .  ~  c . 
_J _ L 


f  . 
J 


if  {fg^hy  <  o 

for  one  or  more  j 


e  =  °°,  if  {f„}T{v. }  -  f .  >  o  ,  for  all  j 
C  B  j  j  — 

T 

Therefore,  by  choosing  £  _<  £  ,  minimizing  z+  =  {c  }  {x}  also  minimizes 
T  ^ 

2  =  {c}  {x}  and,  at  the  same  time,  maximizes  z*  =  E { f }  {x}. 
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APPENDIX  III. 


NOTATION 


A 

[A] 


[b] 

{cb>T 


G 


n 


I 


L 

i 


Mi 

ML 


Mj. 


The  following  symbols  are  used  in  this  paper: 

=  differential  operator; 

=  constraint  matrix; 

=  unknown  coefficients; 

=  columns  of  [a]  matrix; 

=  boundary  differential  expression; 

=  basis  matrix; 

=  row  price  vector; 

=  row  vector  of  the  prices  of  basis  variables 
=  Young's  modulus  of  elasticity; 

=  prescribed  function; 

=  some  specific,  but  arbitrary  row  vector; 

=  shear  modulus ; 

=  nth  order  Hermitian  shape  function; 

=  moment  of  inertia  of  cross  section; 

=  stiffness  of  elastic  foundation; 

=  length  of  a  beam; 

=  length  of  a  beam  element; 

=  moment  at  node  i; 

*  plastic  moment  of  a  beam; 

*  torque ; 

=■  uniform  load  intensity; 


R 

r 

u 

u 


V 

1 


X 

y 


“1 

z  . 

3 

<M*> 


4> 


\ 


v 

a 

V 


residual ; 

maximum  of  the  absolute  value  of  residual 
exact  solution; 

approximate  solution  (trial  function); 

shear  at  node  i  ; 

vector  of  space  variables; 

yield  function; 

value  of  y  at  node  j; 

[B'V*.}; 

9  y 

value  of  -B-tr  at  node  j; 
objective  function; 
secondary  objective  function; 

{cb}T 

basis  of  a  trial  function; 

Prandtl's  stress  function; 

value  of  £  at  node  (i,j)  (Refer  to  Fig.  3) 

3<t> 

value  of  at  node  (i,j); 

J  '-7 

load  factor; 

Poisson's  ratio;  and 

angle  of  twist  per  unit  length. 
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TABLE  1  Results  for  a  Beam  on  Elastic  Foundation  Subject  to  Uniform  Load 
by  using  an  Ordinary  Polynomial 


TABLE  2  Results  for  Displacement  for  Beam  on  Elastic  Foundation  Subject  to  Uniform 
Load  by  using  Hermit  inn  Shape  Function 


(a)  Beam  with  two  elements 

0.305  m;  1  lb/ft  =  14.59  N/m;  1  lb/ft2 


El 
1 

y.n 

(b)  Beam  element 

47.837  N/m2;  1  lb-ft2  =  0.414  N-m2) 


Fig.  1 


Beam  on  Elastic  Foundation 


Fig. 
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25.4  nun;  1  lb  =  4.45  N;  1  lb/in  =  0.175  N/mm) 


Fig.  3  Pure  Tension 


(a)  Original  structure 


(b)  Structure  for  the  first  analysis 


(c)  Structure  for  the  second  analysis 


(d)  Collapse  mechanism 


1  lb-in~  =  2870.692  N-mm“;  1  lb-in  -  113.03 


Limit  Analysis  of  a  Fixed-fixed  Beam 
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